Supersolid phases in a realistic three-dimensional spin model 
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Supersolid phases, in which a superfluid component coexists with conventional crystalline long 
range order, have recently attracted a great deal of attention in the context of both solid helium and 
quantum spin systems. Motivated by recent experiments on 2iif-AgNi02, we study the magnetic 
phase diagram of a realistic three-dimensional spin model with single-ion anisotropy and compet- 
ing interactions on a layered triangular lattice, using classical Monte Carlo simulation techniques, 
complemented by spin- wave calculations. For parameters relevant to experiment, we find a cascade 
of different phases as a function of magnetic field, including three phases which are supersolids 
in the sense of Liu and Fisher. One of these phases is continuously connected with the collinear 
ground state of AgNiCb, and is accessible at relatively low values of magnetic field. The nature of 
this low- field transition, and the possibility of observing this new supersolid phase in AgNiCb, are 
discussed. 
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Solids and liquids are very different. Placed under 
stress, a liquid will flow, while a solid resists deforma- 
tion. The idea of a supersolid, a state which combines 
the properties of a solid with those of a perfect, non- 
dissipative superfluid, seems therefore to fly in the face 
of common sense. None the less, the proposal that a su- 
persolid might occur through the Bose-Einstein conden- 
sation of vacancies in a quantum crystal [lj, was propelled 
to the centre of debate by recent experiments on 4 He[l]. 

A radically different approach to supersolids was ini- 
tiated by Liu and Fisher [3j, who realised that quantum 
magnets could support states which break the transla- 
tional symmetry of the lattice (and are therefore solids) 
while simultaneously breaking spin-rotational symmetry 
within a plane, a form of order analogous to a superfluid. 
It is now well established that models of two-dimensional 
frustrated magnets with easy-axis anisotropy can support 
such supersolid states jj. Moreover, since the states of 
a spin- 1/2 quantum magnet are in one-to-one correspon- 
dence with hard-core bosons, these supersolids might also 
be realised using cold atoms on optical lattices. Nonethe- 
less, candidates for supersolid states among real, three- 
dimensional magnets remain scarce. An interesting sys- 
tem in this context is the triangular easy- axis magnet, 
2i^-AgNi0 2 i. 

AgNi02 is a very unusual material, built of stacked, 
two-dimensional nickel-oxygen planes, held together by 
silver ions. It combines metallicity and magnetism, with 
the magnetic ions in each plane forming a perfect trian- 
gular lattice, nested within a honeycomb network of con- 
ducting sites [5[. In the absence of magnetic field AgNi02 
supports a stripe-like collinear antiferromagnetic ground 
state, illustrated in Fig.QJa). Recently, AgNi02 has been 
shown to undergo a complicated set of phase transitions 
as a function of magnetic field [6]. Of particular interest 
is the transition out of the collinear ground state at low 
temperatures. 
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FIG. 1: (Color online) a) Low-field collinear "stripe" ground state, 
with spins aligned along the magnetic easy axis (z-axis). b) Re- 
lated supersolid phase for magnetic field parallel to the easy axis : 
"down" spins cant into plane perpendicular to field, while "up" 
spin remain aligned with the field, c) First-neighbour Ji, second- 
neighbour J2, and interlayer interactions J± for a stacked a trian- 
gular lattice. 



In applied magnetic field, collinear antiferromagnets 
with easy-axis anisotropy typically undergo a first-order 
"spin-flop" transition into a canted state, at a critical 
field which is broadly independent of temperature. How- 
ever the low-field transition in AgNi02 is accompanied 
by a relatively broad feature in specific heat, does not 
exhibit marked hysteresis, and occurs at progressively 
higher fields as temperature increases. None of these fea- 
tures resemble a typical spin-flop transition, and together 
they raise the question of whether a novel type of mag- 
netic order is realised in AgNi02 under field. 

In this Letter we explore the different phases that occur 
as a function of magnetic field in a simple effective spin 
model already shown to provide excellent fits to inelastic 
neutron scattering spectra for AgNiC^Q- We show that 
the collinear ground state of this model does not undergo 
a conventional spin- flop, but rather a Bose-Einstein con- 
densation of magnetic excitations which converts it into 
it a state that is a supersolid in the sense of Liu and 
Fisher. We also identify two magnetization plateaux, and 
two further supersolid phases at high field. 
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The model we consider is the Heisenberg model on a 
layered triangular lattice, with competing antiferromag- 
netic first- and second- neighbour interactions J\ and J2, 
single-ion anisotropy D and inter-layer coupling J± 



^E,(Sf) 2 -^E,sf. (1) 



(cf. Fig. ffl^c)). For concreteness, we set J± = 1, J2 = 
0.15, J± = —0.15 and D = 0.5, measuring magnetic field 
h and temperature T in units of J\. These are ratios 
of parameters comparable to those used to fit inelastic 
neutron scattering spectra for AgNi02[7]. Like AgNi02, 
in the absence of magnetic field, this model exhibits a 
collinear "stripe-like" magnetic ground state, illustrated 
in Fig.[]Ja). The stripes have three possible orientations, 
and so break a Z3 rotational symmetry of the lattice. The 
collinear stripe state also breaks translational symmetry 
in the direction perpendicular to the stripes. But in the 
presence of a magnetic easy axis, it does not break spin 
rotation symmetry. 

This collinear stripe state supports two branches of 
spin- wave excitations, which are degenerate in the ab- 
sence of magnetic field. Both branches are gapped. 
However for parameters relevant to AgNiC>2, the mini- 
mum of the spin wave dispersion does not occur at the 
magnetic ordering vector M, as would be expected, but 
rather at points M r related by the broken Z3 rotational 
symmetry [7j. 
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FIG. 2: (Color online) (a) Linear spin-wave dispersion of the 
collinear stripe phase of Eq. [T] for h=0 (dashed line) and h=2DS 
(solid lines) in the k z =0 plane, for S = 1, Ji=l, J2=0.15, Jj_=-0.15 
and D=0.5. (b) First Brillouin zone for a triangular lattice, show- 
ing the ordering vector M, and related symmetry points M' . The 
magnetic Brillouin zones for the collinear stripe phase and associ- 
ated supersolid are shown by a blue rectangle and a red hexagon, 
respectively. 

Applying a magnetic field parallel to the easy axis 
lifts the degeneracy of the two spin wave branches, and 
reduces the gap at M' . Within linear spin wave the- 
ory, neglecting dispersion in the out-of-plane direction 
and expanding about a stripe state with ordering vector 
M=(0,2tt/ v / 3), we find 



e±(k) =45 



(ji cos 2 (k x /2) + J 2 cos 2 (V3k y /2) + D/ 2 )' 



(cos(VSk y /2) (Ji cos(k x /2) + J 2 cos(3k x /2))y 



±h. 



and the spin gap at M' =(±7r, ir/y/S) closes completely 
at a critical field h=2DS. The resulting dispersion is 
shown in Fig. [2] (a). 

As in the celebrated example of TICUCI3, the closing 
of this spin gap leads to Bose-Einstein condensation of 
spin- wave excitations (magnons) [i| . This Bose-Einstein 
condensate breaks a 17(1) spin-rotation symmetry in the 
S x -S y plane, and so has superfTuid character. Since the 
resulting state inherits the broken Z2 translational and 
Z3 rotational symmetries of the collinear ground state, 
it is a supersolid. This quantum phase transition can 
also be understood at a mean-field level — instead of 
undergoing a spin- flop, the "down" spins cant, while the 
"up" spins remain aligned with the field. The nature of 
the new magnetic supersolid is illustrated in Fig. QJb). 

These arguments establish the possibility of a super- 
solid state in AgNi02, but tell us nothing about its 
thermodynamic properties. If the low-field transition in 
AgNi02 is into a supersolid, why does the critical field 
increase with increasing temperature ? What might the 
experimental signatures of this new phase be ? What 
other states might occur at higher magnetic field, and 
how do they evolve with temperature ? 

In order to address these questions, we have performed 
classical Monte Carlo simulations of Eq. [TJ The combi- 
nation of large magnetic anisotropy and competing in- 
teractions means that simulations based on a Metropolis 
algorithm with a simple local update suffer severe freez- 
ing. To overcome these problems we employed a par- 
allel tempering Monte Carlo scheme [Io|, combined with 
successive over-relaxation sweeps [Uj. Simulations of 48- 
128 replicas were performed for rhombohedral clusters 
of 3Lx3LxL=9L 3 spins, where L=4,6,8,10 counts the 
number of triangular lattice planes. Periodic boundary 
conditions were imposed. Typical simulations involved 
4xl0 6 steps, half of which were discarded for thermal- 
ization. Each step consisted of one local-update sweep 
of the lattice, followed by two over-relaxation sweeps, 
with replicas at different temperatures exchanged every 
10 steps. We set |S| = S = 1 throughout. 

The results of these simulations are summarised in 
Fig. 03 For the parameters used, we find a total of six 
distinct phases as a function of increasing field : i) a 
collinear "stripe" ground state with a 2-site unit cell; ii) 
a supersolid phase with a 4-site unit cell, continuously 
connected with (i); iii) a collinear one-third magnetiza- 
tion plateau state with a 3-site a/3 x a/3 unit cell; iv) 
a second supersolid, formed by a 2:1:1 canting of spins 
within the same 4-site unit cell as (ii); v) a collinear half- 
magnetization plateau state, with the same 4-site unit 
cell as (ii); and vi) a third supersolid, with the same 
4-site unit cell as (ii), formed by a 3:1 canting of spins 
approaching saturation. Phase transitions were identified 
using peaks in the relevant order-parameter susceptibili- 
ties. These transitions are generically first order, except 
between collinear and supersolid phases with the same 
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FIG. 3: (Color online) Magnetic phase diagram obtained from 
classical Monte Carlo simulation of Eq. ^for a cluster of 24x24x8 
spins with J±=l, J2=0.15, Jj_=-0.15, D=0.5. Temperature T and 
magnetic field h are measured in units of J\ . Phase boundaries are 
determined from peaks in corresponding order parameter suscepti- 
bilities. Phase transitions are first order, except where shown with 
a dashed line. A dotted black line shows the cut at h=1.25 used in 
Fig. HJa)-(d). Inset shows the range of parameters for which a su- 
persolid arises as the first instability of the stripe phase in magnetic 
field, as determined by mean field calculations for J^=-0.15. 



unit cell, where transitions are found to be continuous. 
All of these phases can also be found in mean field the- 
ory at T = 0, and transitions between them are shown by 
open squares, diamonds or circles on the h-axis of Fig. [3l 

This phase diagram shows some intriguing similarities 
with experimental work on AgNiC^Q. In particular, 
the topology of the low-field phases is correctly repro- 
duced, with the low-field supersolid phase contained en- 
tirely within the envelope of the collinear stripe phase. 
The phase transition between these two phases is con- 
tinuous, and the critical field increases with increasing 
temperature [l5| . For this reason we now concentrate on 
the low-field properties of the model, leaving the rich 
physics at higher field for discussion elsewhere. We note, 
however, that the one-third magnetization plateau (hi) 
is well known from studies of easy-axis triangular lattice 
antiferromagnets[12|, and that states analogous to the 
half-magnetization plateau (v) and high field supersolids 
(iv) and (vi) also occur in models of Cr spinels [13j|. 

As in some previously studied models 0, Isl? two finite- 
temperature phase transitions separate the low-field su- 
persolid phase from the paramagnet. The first of these 
is a first-order transition into the collinear stripe state 
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FIG. 4: (Color online) a) U(l) order parameter showing onset 
of supersolid phase for /z=1.25, T^O.l (inset : crossing of associ- 
ated Binder cummulants). b) related order parameter susceptibil- 
ity Xu(i)- c ) ^3 order parameter showing onset of collinear stripe 
phase for h=1.25, T w 0.42 (inset : Binder cumulants for energy, 
showing a dip indicative of a bimodal distribution), d) related order 
parameter susceptibility x^ 3 - Results are taken from simulations 
of clusters of 3Lx3LxL spins, with L=4, 6, 8, 10, for parameters 
identical to Fig. [3] e) finite-size scaling of order parameter suscep- 
tibility at transition into supersolid for h=1.5, T«0.2. f) graphical 
representation of U(l) order parameter as a vector in the S x -S y 
plane. 



at a temperature T « 0.42. The second is a continuous 
transition at a critical temperature which varies approxi- 
mately linearly with magnetic field from T = (h = 1) to 
T « 0.3 (h ~ 2). Both translational and rotational lat- 
tice symmetries are broken at the upper transition. To 
study this it is convenient to introduce a two-component 
order parameter based on an irreducible representation 
of the Cs = Z3 rotation group, which measures the ori- 
entation of the "stripes" in the plane 
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where 5\ = (1,0) and 82 = (1/2, \/3/2) are the primitive 
vectors of the triangular lattice. Figure HJc) and (d) show 
the behaviour of this order parameter, Binder cumu- 
lants for energy, and related susceptibility for h = 1.25, 
T w 0.42. The transition is clearly first order; we have 



checked explicitly that the Z2 symmetry associated with 
translations is broken at the same temperature. 
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FIG. 5: (Color online) Magnetic torque r = m x h as a function of 
temperature, for a magnetic field of magnitude 1.25 at an angle of 
5° to the easy axis (natural units). Torque changes sign abruptly at 
the first order transition from paramagnet to collinear stripe phase 
for T « 0.42. A change in slope for T « 0.1 signals the continuous 
transition from stripe phase to magnetic supersolid. Insets show 
the heat capacity anomalies at each of these transitions. 



Spin rotation symmetry in the S x -S y plane is broken 
at the lower phase transition into the supersolid state. 
This can be measured by constructing a (7(1) order pa- 
rameter which measures the difference between S x and 
S y components of the canted "down" spins, as illustrated 
in Fig. 0(f). Fig. 2(a) and (b) show the behaviour of this 
order parameter, its Binder cumulants, and related sus- 
ceptibility for h = 1.25, T & 0.1. The phase transition 
remains continuous at finite temperature, with Binder 
cumulants for different system size crossing at a single 
temperature. A good collapse of susceptibility data is 
obtained using susceptibility and correlation length ex- 
ponents 7 = 1.32 and v = 0.67 for the 3D XY universality 
class, as shown in Fig. E(e). 

While the relative extent of each phase and details of 
critical fields and temperatures are different, AgNiC>2 ex- 
hibits a similar double transition on cooling; a transition 
from paramagnet to a collinear stripe phase at T^20K 
accompanied by a sharp feature in specific heat and then, 
for fields greater than 13.5T, a continuous or very weakly 
first-order transition from the collinear stripe phase into 
an unknown low temperature magnetic state. This tran- 
sition occurs at higher magnetic field for higher temper- 
atures, suggesting that the stripe phase has higher en- 
tropy than the competing high-field phase, as found in 
our simulations. Is the high field phase in AgNiC>2 then 
a supersolid ? 

Direct confirmation of the magnetic order in AgNi02 
for fields greater than 13.5T by elastic neutron scattering 
is feasible, but challenging, since no large single crystals 
are presently available. None the less, it should be possi- 
ble to observe the closing of the spin gap on entry to the 
supersolid phase in inelastic neutron scattering experi- 
ments on powder samples [7]. Moreover, both transport 



and thermodynamic measurements on small single crys- 
tals clearly resolve magnetic phase transitions in AgNi02 
(sj . We therefore conclude by examining the heat capac- 
ity and magnetization (torque) signatures of the stripe 
and supersolid phases of our model. 

In Fig. [5] we present predictions for magnetic torque, 
r = m x h, and heat capacity Ch, for the same fixed 
value of magnetic field h = 1.25 chosen for the study 
of phase transitions in Fig. |4j Torque changes sign at 
the first-order phase transition from the paramagnet at 
T « 0.4, is strongly temperature dependent in the stripe 
phases for 0.1 < T < 0.4, and is broadly temperature- 
independent in the supersolid phase for T < 0.1 [Hj. For 
the ratios of parameters used, with Ji=1.32meV (cf.Q), 
this translates into a supersolid transition at a field of 
12.5 Tesla, for a temperature of 1.5K. The heat capac- 
ity anomalies at both transitions strongly resemble those 
observed in AgNiC^lfil- 

In summary, we have studied the magnetic phase di- 
agram of a realistic three-dimensional spin model with 
single-ion anisotropy and competing interactions on a 
layered triangular lattice, identifying three phases which 
are magnetic supersolids in the sense of Liu and Fisher 
We find that these supersolids are continuously connected 
with parent collinear phases through the Bose-Einstein 
condensation of magnons. However, quantum fluctuation 
effects may also play an important role at these phase 
transitions, and this remains an interesting problem for 
future study. The model studied was motivated by the 
metallic triangular lattice antiferromagnet 2i^-AgNi02, 
and is known to describe its magnetic excitations in zero 
fieldQ. Since the model does not take itinerant charge 
carriers into account, it cannot pretend to be a com- 
plete theory of AgNi02- Nonetheless, it motivates a 
re-examination of the low field transitions observed in 
AgNi02, where a magnetic supersolid may already have 
been observed Q. 
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